Breakdown of Universality in Quantum Chaotic Transport: 
the Two-Phase Dynamical Fluid Model 



Ph. Jacquod and E.V. Sukhorukov 

Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve 4, Switzerland 

(Dated: November 22, 2003) 

We investigate the transport properties of open quantum chaotic systems in the semiclassical 
limit. We show how the transmission spectrum, the conductance fluctuations, and their correlations 
are influenced by the underlying chaotic classical dynamics, and result from the separation of the 
quantum phase space into a stochastic and a deterministic phase. Consequently, sample-to-sample 
conductance fluctuations lose their universality, while the persistence of a finite stochastic phase 
protects the universality of conductance fluctuations under variation of a quantum parameter. 
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Universal Conductance Fluctuations (UCF) are ar- 
guably one of the most spectacular manifestations of 
quantum coherence in mesoscopic systems |l| . In metallic 
samples, the universality of the conductance fluctuations 
manifests itself in their magnitude, rms(<?) = 0(e 2 /h), 
independently on the sample's shape and size, its average 
conductance or the exact configuration of the underlying 
impurity disorder 0, 0] ■ In ballistic chaotic systems, a 
similar behavior is observed .which is captured by Ran- 
dom Matrix Theory (RMT) [3j . For an open chaotic cav- 
ity connected to two iV-channel leads, and thus having 
an average classical conductance g — N/2 (we consider 
spinless fermions and express g in units of e 2 /h), RMT 
predicts a universal conductance variance cr 2 (g) = 1/8 
for time-reversal and spin rotational symmetric samples. 
At the core of the UCF lies the ergodic hypothesis that 
sample-to-sample fluctuations are equivalent to fluctua- 
tions induced by parametric variations (e.g. changing the 
energy or the magnetic field) within a given sample p) . 

According to the scattering theory, transport proper- 
ties such as the conductance derive from transmission 
eigenvalues, i.e. g = 5Z i=1 Tj B- While coherence ef- 
fects such as the UCF arise due to nontrivial correla- 
tions (TiTj), the knowledge of the probability distribu- 
tion P(T) of transmission eigenvalues is sufficient to cor- 
rectly predict, e.g., the average conductance, or the Fano 
factor F = (T(l - T))/(T> for the shot-noise power [|. 
For a ballistic chaotic cavity, RMT predicts 

Prmt(T) = - 1 (1) 

7T y/T(l-T) 

and thus F = 1/4. For shot- noise, as well as for UCF, 
the correct universal behavior is captured by RMT. 

The validity of RMT is however generically restricted 
by the existence of finite time scales. Spectral fluctu- 
ations are known to deviate from RMT predictions for 
energies larger than the inverse period of the shortest 
periodic orbit for chaotic systems [fj|, or than the in- 
verse time of diffusion through the sample in disordered 
metallic systems Q. Another time scale which is ab- 
sent in RMT is the Ehrenfest time te @> i-e. the time 
it takes for the underlying classical chaotic dynamics 
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FIG. 1: Distribution of transmission eigenvalues for K = 
27.65, td = 25, and M = 2048 (empty circles; te ~ 0; dis- 
tribution calculated over 729 different samples); K = 9.65, 
td = 5, and M = 1024 (black diamonds; te = 1.5; 729 sam- 
ples); M = 8192 (empty squares; te = 2.8; 16 samples); and 
M — 65536 (black triangles; te ~ 4.1; 1 sample). The solid 
line gives the universal distribution Prmt of Eq. fT| . and the 
dashed line the distribution P a of Eq. with a = 0.39. 
Note that P(T) is symmetric around T = 0.5. Inset: in- 
tegrated probability distribution of transmission eigenvalues 
for the same set of parameters as in the main panel, as well 
as for K = 9.65, t d = 5, and M = 128 (x; t e = 0.16). The 
solid lines are fits obtained from Eq. with a ~ 0.98, 0.81, 
0.6, 0.45 and 0.385 (from bottom to top). 



(with Lyapunov exponent A) to stretch an initial nar- 
row wavepacket, of spatial extension given by the Fermi 
wavelength \p, to the linear system size L. Defining 
M = L/X F , one has t e = A -1 \u[M/(2t d ) 2 } HQ3, with 
Td = M/2N, the dwell time through the cavity (all times 
will be measured in units of the time of flight across the 
cavity). Note, in particular, that the growth of te in 
the semiclassical limit M — > oo is only logarithmic. The 
emergence of a finite t e /td leads to strong deviations 
from the universal RMT behavior, and in particular to 
the suppression of shot noise [lllll2L ll3i . or the proximity 
gap in Andreev billiards 0, Il4l Il5l Il6j| . It has further- 
more been predicted that weak localization vanishes at 
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large te/td [l3- Also, in dirty d-wave superconductors, 
the RMT behavior of the quasiparticle density of states 
[l8l | is restored only below an energy scale set by te 0] • 
The suppression of shot-noise for te/te — > oo is due 
to the disappearance of the stochasticity of quantum me- 
chanical transport, and its re plac ement by the determin- 
ism of classical transport [llL Il2l Ha ] . Wavepackets trav- 
eling on scattering trajectories shorter than te have no 
time to diffract, and are thus either fully transmitted or 
fully reflected. In an open chaotic cavity, different scat- 
tering trajectories have in general different dwell times 
with a distribution p(t) = exp(— i/r^/ru. For finite 
< t e /t d <C oo, this suggests that transport is mediated 
by a two-phase dynamical fluid, consisting of a stochastic 
phase of relative volume a ~ / p(t)dt = exp(— te/te), 
and a deterministic phase of relative volume 1 — a. Fol- 
lowing this purely classical argument, first expressed in 
Ref. [l2j, one expects that the distribution of transmis- 
sion eigenvalues is given by 



P a (T) = aPRMT(T) 



1-a 



[6(T) + S(l T)] . (2) 



This will be confirmed below. We will further show that, 
quite surprisingly, the two-phase dynamical fluid assump- 
tion also correctly describes the behavior of mesoscopic 
coherent effects, and in particular that it explains the 
breakdown of universality of the conductance fluctua- 
tions when te becomes comparable to the ergodic time 

TO- 

We first summarize our main results, (i) We give 
full confirmation of Eq. J5J by calculating the inte- 
grated distribution of transmission eigenvalues I(T) = 
J T P(T')dT'. We find that it is very well fitted by (see 
the inset to Fig.l) 
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from which we extract a ~ exp(— te/te)- (ii) The con- 
ductance fluctuations stay at their universal value, inde- 
pendently on te/td, under variation of the energy in a 
given sample. This follows from the survival of a large 
number of stochastic channels - even though their rela- 
tive measure a — » - which preserves the universality of 
the conductance fluctuations, (iii) A completely differ- 
ent situation arises when one considers sample-to-sample 
fluctuations. In this case, one has <J 2 {g) oc (M/M c ) 2 for 
M > M c . The scaling parameter M c = t d exp(A) is de- 
termined by the quantum mechanical resolution of clas- 
sical phase space structures corresponding to the largest 
cluster of fully transmitted or reflected neighboring tra- 
jectories (see Ref. 12]). (iv) The energy conductance 
correlator always decays on the universal scale of the 
Thouless energy, £ £ oc 1/td, independently on te- The 
results (ii) and (iii) show that the ergodic hypothesis 
breaks down as t e /t increases. Accordingly, (iv) is 
somewhat surprising, but will be understood below via 
a semiclassical argument. All our results and arguments 



fully confirm the two-phase dynamical fluid model. We 
note that our conclusion (iii) is in agreement with the 
very recent finding a 2 (q) oc M 2 obtained by Tworzydlo, 
Tajic and Beenakker |20|. However, (ii) is in complete 
opposition with their prediction that deviations from 
cr 2 (g) = 1/8 should occur upon variation of the energy 
for large te/te- Points (i) and (iv) are addressed here 
for the first time. 

We consider open systems with fully developed chaotic 
dynamics, for which td ^ 1- Because te grows logarith- 
mically with M, and since we want to investigate the 
regime te/te }t 1) wc model the electron dynamics by 
the kicked rotator map, which reproduces most of the 
phenomenology of low-dimensional noninteracting elec- 
tronic physics |22|. The classical kicked rotator map is 
given by 



x = x + p 

p = p + K sin(x), 



(4) 



with K the (dimensionless) kicking strength. It drives the 
dynamics from fully integrable (K = 0) to fully chaotic 
[K > 7, with Lyapunov exponent A ps ln(K/2)]. We 
consider a toroidal classical phase space x,p S [0, 2ir], and 
open the system by defining contacts to ballistic leads via 
two absorbing phase space strips [xl — Sx, xl + Sx] and 
[xr — Sx, xr+5x], each of them with a width 2Sx = tt/te- 
Quantizing the map amounts to a discretization of, 
say, the real space coordinates as x m = 27rm/M, m = 
1,...M. A quantum representation of the map 01 is 
provided by the unitary M x M Floquet operator U [21| , 
which gives the time evolution for one iteration of the 
map. For our specific choice of the kicked rotator, the 
Floquet operator has matrix elements 
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fyf—l /2 e —(iM K I ^47r)[cos(2-irm / " M)+cos(2nm' / 'M)\ 
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2TTil(m-m')/M e -(ni/2M)l 2 



(5) 



The spectrum exp(ie Q ) of U defines a discrete set of M 
quasienergies e a € [0, 2ir) with an average level spacing 
S = 2-n/M. 

Much in the same way as in the Hamiltonian case [2^] , 
a quasienergy-dependent 2N x 2N scattering matrix can 
be determined from the Floquet operator U as 24] 

S{e) = P[exp(-ie) - U(l - P T P)]- 1 UP T , (6) 

using a 2A^ x M projection matrix P which describes the 
coupling to the leads. Its matrix elements are given by 



1 if n = m G {m 
otherwise. 



(L) 



}U{mf ] } 



(7) 



An ensemble of samples with the same microscopic prop- 
erties can be defined by varying the position {m[ L ' R ^}, 
i = 1, ... N of the contacts to the left and right leads for 
fixed T D = M/2N and K. 
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FIG. 2: Variance u 2 (g) of the conductance vs. M/M c , for 
microscopic parameters K G [9.65,27.65], td G [5,25], and 
M G [128, 16384]. The scaling parameter M c t d exp(A) varies 
by a factor 70. The solid and dashed lines indicate the clas- 
sical, sample-to-sample behavior oc M 2 , and the universal 
behavior u 2 {g) = 1/8 respectively. Inset: unsealed data for 
K — 9.65 and td = 5 (circles), 7 (squares), 10 (diamonds), 
15 (upward triangles) and 25 (downward triangles). 

As usual, the scattering matrix can be written in a four 
block form in term of N x N transmission and reflection 
matrices as 

S = (J'/)" (8) 

The spectrum of transmission probabilities is given by 
the N eigenvalues Ti of T = ffi , from which the dimcn- 
sionless conductance is obtained, via the Landauer for- 
mula g = Ti Our numerical procedure follows the 
description given in Ref. |l3[ . 

We plot in Fig.l various distributions P(T) of trans- 
mission eigenvalues. First, it is seen that our model cor- 
rectly reproduces the RMT distribution of Eq. in the 
limit te/tjj -C 1. The distribution undergoes strong 
modifications, however, as te/tjj increases. In particu- 
lar, more and more weight is accumulated at T — and 1. 
The behavior exhibited by P(T) for finite te/tq seems 
very similar to that predicted by Eq. @. To confirm 
this, we calculate the integrated probability distribution 
I(T), which presents the advantage of being independent 
on the size of histogram bins. Results are shown in the 
inset to Fig. 1. The fitting curves clearly confirm the va- 
lidity of Eq. (J2J) . The extracted parameter a is found to 
obey a « exp(— (1 + te)/td), for te > 0. We attribute 
the factor 1 + te in the exponential (and not Te) to the 
discrete nature of the dynamics in our model. 

For te — > 0, one is in the UCF regime, where the con- 
ductance fluctuates equivalently from sample to sample 
or as e is varied within a given sample. This is no longer 
the case, however, once te becomes finite, as is shown in 
Fig. 2. While a 2 (g) = 1/8 seems to be preserved when 
e is varied for a given sample, one gets an enormous in- 
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FIG. 3: (a) Phase space cross section of the left lead for 
K = 9.65, and td — 5. Black dots indicate transmitted, 
white areas reflected classical trajectories respectively. We 
used 25000 initial conditions, and only trajectories exiting the 
system after less than 5 iterations of the classical map @ have 
been kept, (b) Correlation length £ e extracted from the con- 
ductance correlator as F(£ s ) = 0.8, for K = 9.65 and td = 5 
vs. M (black circles) and for K = 9.65 and M = 2048 vs. 
td (white squares). The dashed line indicates the expected 
~ 1/td behavior (see text), (c) Decay of the conductance 
correlator F(e) vs. quasienergy for K = 9.65, M = 2048, and 
td = 5, 10, 15, 20, and 40. 

crease <J 2 (g) oc M 2 from sample to sample. This behav- 
ior derives from the underlying classical dynamics, and 
can be understood on the basis of a two-phase dynamical 
fluid, as we now proceed to explain. 

While the classical dynamics considered here is fully 
chaotic, finite-sized phase space structures emerge due 
to the opening of the cavity. Following Ref. these 
structures can be visualized by marking which trajecto- 
ries originating from, say, the left lead, ends up being 
transmitted to the right lead. Such a picture is shown in 
Fig. 3(a). It is seen that classical trajectories are trans- 
mitted by a series of bands exiting the cavity at times 
tj, and covering a phase space area Aj — t^ 2 exp(— Xtj). 
The shape, position and precise volume of these bands is 
sample-to-sample dependent. In the semiclassical limit 
fi c s — 2-k/M — * 0, the effective Planck scale Ti c g resolves 
the j th band as soon as h e s < Aj, or M > 2ttte exp(Atj-). 
Once the largest classical band is resolved, a 2 (g) starts 
to be dominated by the band fluctuations. Each resolved 
band carrying a growing number oc M of fully transmit- 
ted (or reflected) quantum mechanical modes, one ex- 
pects a variance <r 2 (g) oc (M/M c ) 2 for sufficiently large 
M > M C) with a scaling parameter M c = InTp, exp(A) 
determined by the largest band, exiting the system at the 
ergodic time tq « 1. As shown on Fig. 2, this is precisely 
what happens pEj . Deviations from the universal behav- 
ior emerge for M w M c , equivalently when te/tq ~ 1, 
i.e. much earlier than the suppression of shot noise (see 
also Ref. [U). 
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We face a completely different situation when varying 
e within a given sample. In a ballistic system like ours, 
such a change does not modify the classical trajectories, 
and thus alters only the action phase of each contribution 
to the semiclassical Green function. Within this picture, 
the conductance fluctuates only due to long, diffracting 
orbits with t > te 26]. These long orbits build up the 
stochastic phase. Their subset can be viewed as corre- 
sponding to an effective stochastic cavity with contacts to 
leads with N e g — aN channels. Fixing the microscopic 
parameters t d and A, one has N e g ~ Af 1_1 / ATB 1. 
This means that, despite the prefactor a, N c g is always 
large enough to guarantee that transport occurs semiclas- 
sically, and therefore, one stays always in a regime with 
universal value u 2 (g) = 1/8 3]. A first confirmation of 
this argument is provided by the numerical data shown 
in Fig. 2, which indicate a constant behavior of o~ 2 (g), in- 
dependently on te/td- To further check this argument, 
we finally consider the conductance correlator 

F(e) = a- 2 (g) (Sg(s )Sg(e + e)) . (9) 

As said above, only the phase accumulated after diffrac- 



tion (for t > te) contributes to conductance fluctuations 
|2bT |. and since the subset of diffractive trajectories have 
an average dwell time given by te + Trj, they accumulate 
a relevant relative phase oc s Tp. One therefore expects 
a decay of F(e) over the Thouless scale as in the univer- 
sal regime [ij, £ e oc 1/td, independently on te- This is 
confirmed by the data shown in Fig. 3(b) and (c). 

Our results thus show that the separation of the de- 
terministic and stochastic phases is complete. Beyond a 
simple explanation of the suppression of shot-noise with 
a = exp(— te/td) via Eq.Q, the phase separation cor- 
rectly accounts for the behavior of the conductance vari- 
ance and correlators in open quantum chaotic systems 
in the semiclassical limit. Further investigations along 
the lines initiated here should focus on other effects of 
mesoscopic coherence, such as the weak-localization cor- 
rections, and in particular the magnetoresistance. 
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National Science Foundation. 



[1] B.L. Altshuler, JETP Letters 41, 648 (1985); P.A. Lee 
and A.D. Stone, Phys. Rev. Lett. 55, 1622 (1985). 

[2] A. Benoit, CP. Umbach, R.B. Laibowitz, and R.A. 
Webb, Phys. Rev. Lett. 58, 2343 (1987); W.J. Skocpol, 
P.M. Mankiewich, R.E. Howard, L.D. Jackel, D.M. Ten- 
nant, and A.D. Stone, Phys. Rev. Lett. 58, 2347 (1987). 

[3] H.U. Baranger and P.A. Mello, Phys. Rev. Lett. 73, 
142 (1994); R.A. Jalabert, J.-L. Pichard, and C.W.J. 
Beenakker, Europhys. Lett. 27, 255 (1994). 

[4] M. Biittiker, IBM J. Res. Dev. 32, 317 (1988). 

[5] Ya.M. Blanter and M. Biittiker, Phys. Rep. 336, 1 
(2000). 

[6] M.V. Berry, Proc. R. Soc. Lond. A 400, 229 (1985). 
[7] B.L. Altshuler, and B.I. Shklovskii, Sov. Phys. JETP 64, 
127 (1986). 

[8] G.P. Berman and G.M. Zaslavsky, Physica A 91, 450 

(1978) ; M.V. Berry and N.L. Balasz, J. Phys. A 12, 625 

(1979) . 

[9] M.G. Vavilov and A.I. Larkin, Phys. Rev. B 67, 115335 
(2003). 

[10] The deviation from the closed sample result te = 

A -1 ln(M) is due to partial reflection at the boundary 

between the cavity and the contacts to the leads [jj. 
[11] O. Agam, I. Aleiner and A. Larkin, Phys. Rev. Lett. 85, 

3153 (2000); S. Oberholzer, E.V. Sukhorukov, and C. 

Schonenberger, Nature 415, 765 (2002). 
[12] P.G. Silvestrov, M.C. Goorden, and C.W.J. Beenakker, 

Phys. Rev. B 67, 241301 (2003). 
[13] J. Tworzydlo, A. Tajic, H. Schomems, and C.W.J. 

Beenakker, Phys. Rev. B 68, 115313 (2003). 
[14] A. Lodder and Yu.V. Nazarov, Phys. Rev. B 58, 5783 

(1998); D. Taras-Semchuk and A. Altland, Phys. Rev. B 

64, 014512 (2001). 



[15] Ph. Jacquod, H. Schomerus, and C.W.J. Beenakker, 
Phys. Rev. Lett. 90, 207004 (2003); M.C. Goorden, Ph. 
Jacquod, and C.W.J. Beenakker, Phys. Rev. B 68 220501 
(R) (2003). 

[16] A. Kormanyos, Z. Kaufmann, C.J. Lambert, and J. 

Cserti, cond-mat/0309306 
[17] I.L. Aleiner and A.l. Larkin, Phys. Rev. B 54, 14423 

(1996); i. Adagideli, Phys. Rev. B 68, 233308 (2003). 
[18] A. Altland, B.D. Simons, and MR. Zirnbauer, Phys. 

Rep. 359, 283 (2002). 
[19] i. Adagideli and Ph. Jacquod, Phys. Rev. B 69, 

020503(R) (2004). 
[20] J. Tworzydlo, A. Tajic, and C.W.J. Beenakker. 

cond-mat/0311283 
[21] L.E. Reichl, The Transition to Chaos, Springer, New 

York (1992). 

[22] S. Fishman, D. R. Grempel, and R. E. Prange, Phys. Rev. 
Lett. 49, 509 (1984); A. Altland and M. R. Zirnbauer, 
Phys. Rev. Lett. 77, 4536 (1996). 

[23] S. Iida, H.A. Weidenmiiller, and J. A. Zuk, Phys. Rev. 
Lett. 64, 583 (1990). 

[24] Y.V. Fyodorov and H.-J.Sommers, JETP Letters 72, 
422 (2000); A. O ssipov, T. Kottos, and T. Geisel, 
cond-mat/0208378 

[25] The estimate Aj ~ exp( — At,) with = 2 strictly ap- 
plies when the beam of electrons injected into the cavity 
has a similar spread in momentum and real-space coor- 
dinates. For the data shown in Fig. 2, we found a sig- 
nificantly better numerical rescaling with f3 « 1.5 (not 
shown) , which we relate to the anisotropy in the electron 
injection beam. 

[26] N. Argaman, Phys. Rev. B 53, 7035 (1996). 



